Solitons in Trapped Bose-Einstein condensates in one-dimensional optical lattices 
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We use Quantum Monte Carlo simulations to show the presence and study the properties of 
solitons in the one dimensional soft-core bosonic Hubbard model with near neighbor interaction in 
traps. We show that when the half-filled Charge Density Wave (CDW) phase is doped, solitons 
are produced and quasi long range order established. We discuss the implications of these results 
for the presence and robustness of this solitonic phase in Bose-Einstein Condensates (BEC) on one 
dimensional optical lattices in traps and study the associated excitation spectrum. The density 
profile exhibits the coexistence of Mott insulator, CDW, and superfluid regions. 
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Localization phenomena in reduced dimensionality 
take place in a wide range of systems. Shape preserv- 
ing excitations like lattice solitons, which are intrinsi- 
cally localized modes, are caused by an interplay of the 
discreteness of the lattice and non-linearity of the un- 
derlying dynamics. These objects have been the focus 
of intense experimental activity and have been observed 
in optical waveguide arrays in the quasi one dimen- 
sional antiferromagnetic material (C2H5NH3)2CuCl4 
and in quasi-one-dimensional Bose-Einstein condensates 
(BEC) In this last example, a = tt phase flip is 
imposed on a BEC with repulsive interactions (positive 
scattering length), exciting a dark soliton. Topological 
excitations in lattice models of polyacetylene also have a 
long history 0. There, localized regions exist in which a 
transition occurs between two possible configurations of 
long and short bonds. For reviews, see 



Recent effort has focussed on the existence of such 
localized modes in BEC. Dark solitons were exhibited 
in numerical solutions to the Gross-Pitaevskii equa- 
tion 13, while bright solitons were shown to exist 
in BEC with attractive interactions (negative scattering 
length) by solving numerically the non-linear Schrodinger 
equation In addition, bright solitons are known to 
exist for repulsive interactions when the effective mass is 
negative Q and were shown to exist experimentally for 
positive scattering length atomic condensates on trapped 
optical lattices |lOj. Since condensate interactions and 
lattice parameters can be precisely tuned, optical lattices 
offer the possibility to explore systemati cally exotic soli- 
ton phases. Variational calculations [llL have begun 
to explore the dynamics and excitations of such systems. 



In this paper we use Quantum Monte Carlo (QMC) 
simulations to determine the effect of near neighbor (nn) 
repulsive interactions on the ground state phase diagram 
of BEC on Id optical lattices, both with and without 
traps. Our model is the bosonic Hubbard tight binding 



model, 

H = -t^(^ala,+i+al^^a,'j +VT^xfn, (1) 
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The hopping parameter, t, sets the energy scale, rii — 
alai is the number operator, [ai,aj] — Sij are bosonic 
creation and destruction operators. Vr sets the confining 
trap curvature, while the contact and near-neighbor in- 
teractions are given by U and Vi . We use the World Line 
algorithm and the Stochastic Series Expansion (SSE) 
method, which work in the canonical and grand canon- 
ical ensembles, respectively, and simulate both the soft 
core ([/ finite) and hard core cases. 

The phase diagram with Vr = is known at half fill- 
ing ^3 ■ For Vi < 2t the ground state is superfiuid while 
for large U and Vi > 2t off diagonal long range order is 
replaced by an incompressible, insulating charge density 
wave (CDW) phase where sites alternate between high 
and low occupation. Away from half and integer filling 
the system is always superfluid. In the hardcore limit, 
this model can be mapped onto a spinless Id fermionic 
model with near-neighbor repulsion and also onto the 
quantum spin- 1/2 XXZ model. The extended Id fermion 
Hubbard model (with spin) and the classical spin chain 
model are both known to have soliton excitations 2, 1^. 
The main focus of this paper is to demonstrate that soli- 
tons continue to exist despite the possibility of multiple 
occupancy and the presence of a trap. 

Static and dynamic quantities like the density and 
compressibility have already been shown to exhibit un- 
usual features due to the trap, requiring local generaliza- 
tions of these global quantities, [la, [ig We also show that 
even more complex spatial structures can exist when the 
near-neighbor repulsion is nonzero. 

To address the question of solitonic excitations, we 
measure the structure factor at equal imaginary time. 
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where L is the number of lattice sites and < t < (3. Our 
simulations are done at = 10 which is large enough 
to study the ground state. To make contact with the 
excitation spectrum, we use the /-sum rule 



duj ujS{k,uj) = NbEk, (3) 



where S{k,Lu) is the dynamic structure factor (Ni,S{k) 
J dujS{k, to)) and 



L 

(cos k~l) (0| + |0), (4) 

i=l 



The dispersion relation is given by the Feynman result, 
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The dispersion curves shown below are obtained with 
Eq. |31 However, we have verified that we obtain 
the same results by measuring the imaginary-time- 
separated density-density correlation function, perform- 
ing the Laplace transform using the maximum entropy 
algorithm to obtain S{k,uj), and applying Eq. 13 directly. 

We first address the uniform, Vr — 0, system in the 
hardcore limit using the SSE algorithm. Figure ^ shows 
S{k) and ^l{k) for L = 128 sites and two values of Ni, off 
half filling. Note that for small k, ^l{k) for both fillings 
behaves linearly, indicating phonon excitations and su- 
perfluid stability under the Landau criterion. The veloc- 
ity of sound can be extracted from this small k behavior. 
Furthermore, the peak at S{k*), observed at half filling 
for k* — TT, and which indicates the CDW order, does 
not vanish when the system is doped. Instead, its height 
decreases and k* shifts to lower values. Corresponding 
to this quasi-long range order peak in S{k) is a dip in 
ri(fc — k*) giving the soliton excitation energy. This is 
similar to the roton minimum in two and three dimen- 
sions. As the system is doped further, the soliton min- 
imum, and thus quasi-long range order, will disappear. 
The solitons are also evident in the real space boson den- 
sity (not shown) in the form of well-localized regions of 
cross over between the two possible sublattices holding 
high and low density sites. The presence of soliton exci- 
tations in the hardcore system is in agreement with what 
is known for fermionic and classical spin chains. 

The crucial question is whether this behaviour persists 
for the soft core case which is relevant to atomic conden- 
sates on optical lattices. In the soft core case, the contact 
interaction, U, must be large enough to suppress multiple 
occupancy in order to stabilize the CDW phase at half 
filling when Vi is large enough. 13] Such large values of U 
have been achieved experimentally on optical lattices and 
lead to the Mott phase at full filling. In what follows we 
fix [/ = 5t. At Vi = it, the ground state density profile 
and correlation function at half filling exhibit a strong 
CDW pattern. The density-density correlations oscillate 
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FIG. 1: The dispersion relation, Q{k) and the static struc- 
ture factor, S{k), vs k. Peaks in S{k) at incommensurate 
wavevectors result in the soliton minima in Q{k). 




FIG. 2: The density profile and correlation function (aver, 
over 10® QMC sweeps) for the doped soft core system. 



with nearly maximal amplitude, indicating quantum fluc- 
tuations are small, and show little decay with increasing 
separation. Doping this system by removing two bosons 
(Fig. [21 yields pronounced, long-lived, soliton excitations. 
In real space, as seen in Fig. 2, these appear as local 
regions of density alternation modulated by an overall 
'beating' pattern. The beat wavelength (soliton size) is 
given by Ax = 2tt / (n — k*) , where fc* is the position of 
the soliton energy minimum in fl{k). In Fig.|3|we show 
the dispersion ri(fc) for several fillings. S{k) behaves like 
the hardcore case with a peak at k — k* correspond- 
ing to the soliton minimum which moves towards lower 
k* as doping is increased. As in the hard-core case, for 
Nh < 32, where the system is superfluid, fl(k) cx k for 
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FIG. 3: The dispersion relation, fl{k), vs k for several fillings 
for the soft core model, (7 = 5, Vi = 4 at /3 = 10 and L — 64. 
Solitons survive despite multiple occupancy. 



FIG. 5: The excitation energy, fl{k), vs k for the systems 
shown in Fig. 0| For A'^;, — 16 there is a soliton excitation 
at < TT as in the doped uniform system while for A^;, = 22, 
Q{k* = tt) —» indicating the establishment of long range 
CDW order. 




FIG. 4: Density profiles in a trap: For all fillings we see CDW 1^ 
oscillations. For Nt — 22 the central region has long range 

CDW order while Nt = 55 the central region is a MI. The FIG. 6: Evolution of the excitation energy, Q{k), with in- 
A'';, = 16 case exhibits solitonic oscillations. creasing Vi for U — 5t, Nb — 16, Vr = 0.008. The presence 

of solitonic excitations is clear for Vi = 3t and 4t. 



small k showing the stability of the superfluid and the 
absence of a gap. However, for Nt = 32, the system is 
a gapped CDW insulator. This is seen clearly in Fig. |3| 
where r2(fc) goes to a finite value as k ^ 0. There is also 
no soliton feature at intermediate k*. For Ni, = 32 the 
dispersion Q{k) has minima only at fc = 0,7r. 

Placing the system in a trap destroys translational in- 
variance. Nonetheless, we shall now show that robust 
CDW and solitonic excitations are observed. In Fig. ^ 
the local density profiles in a trap, Vr = 0.008, are given 
for three fillings. For iVj, = 16, solitonic oscillations are 
again evident (see also Fig.|3J) as local CDW correlations 
modulated by a beat envelope. For Nf, = 22 long range 



CDW order dominates although some residual solitonic 
excitations remain near the edges. For Nj, — 55 one sees 
a remarkable co-existence of several phases: CDW to- 
wards the edges, followed by superfluid (no CDW oscil- 
lation and compressible) and then a central incompress- 
ible Mott insulator (MI). The density fluctuations in the 
two CDW regions are decoupled by the intervening MI. 
Such striking density oscillations have been observed in 
non-neutral plasmas 0, which, due to their electric 
charge, have long range repulsive interactions. 

Figure shows 0(fc) for the same fillings as in Fig. 0] 
For Ni, = 16 there is a clear solitonic excitation of the 
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FIG. 7: Evolution of the excitation energy, with in- 

creasing Vi for f/ = 5, iVt = 22, I/t = 0.008. The gradual 
establishment of long range CDW order is clear. 

type seen for the uniform system, fc* < tt. For the higher 
fiUing, iVb = 55, the excitation has shifted towards fc* = 
TT but ri(A:) remains relatively high, indicating that this 
is not true long range order (as is of course clear from 
the density profile). For TV^ = 22, on the other hand, 
we see that ^(k* — tt) is close to zero indicating the 
establishment of long range CDW order. Furthermore, 
there is no evidence of a gap in r2(fc) as is seen in the 
uniform system at half filling (see F ig. El ) : The system as 
a whole is always compressible 0, Il6l|. 

Finally, the evolution of the dispersion relation with 
increasing near-neighbor repulsion is shown for N^, = 16 
in Fig. El and — 22 in Fig. [3 In the former case, a 
soliton minimum develops, while in the latter case CDW 
formation takes place instead. 

A further interesting feature of Fig. |B1 with its soli- 
tonic excitations, is the universal crossing of the dis- 
persion curves for different interaction strengths. On 
the other hand, the crossing in Fig. [71 where CDW 
order dominates, is not universal. A similar well de- 



fined crossing point in the specific heat has been seen 
both experimen tally in '^He ,20j and in fermion Hubbard 
mo dels .1211 123 . l23l | We believe a similar reasoning for 
the existence of crossing applies here. The integral of the 
structure factor over all momentum is constrained by the 
density. Thus if S'(fc) increases with Y\ for some momenta 
(for example at A: = tt as CDW correlations build up), 
there must be a corresponding decrease in S(}t) for other 
momenta. This implies a similar behavior in the dis- 
persion relation 0(fc) and hence suggests that dispersion 
curves for different interaction strengths should cross. As 
discussed in |22| the universality of the crossing in the 
specific heat case is a second, and more subtle issue. 

In conclusion, we have demonstrated that soliton sig- 
natures, which are to be expected in the d = \ hard-core 
boson system owing to its close connection with lattice 
fermion and spin models, are still very robust when the 
bosons become soft-core and when they are placed in a 
confining potential. Rapid progress in the creation of 
near neighbor repulsion Y\ in optical trap systems sug- 
gests that it will soon be possible to look for these solitons 
experimentally. 

We have also found that trapped bosons with near 
neighbor repulsion can exhibit a remarkably rich density 
profile in which a Mott insulator at commensurate fill- 
ing occupied the trap center, followed by a superfiuid re- 
gion and then a CDW region where the density is locally 
pinned at ^, with a second, and final superfiuid region 
at the end of the occupied sites. The local compress- 
ibility [isl 0| also exhibits some unusual features. We 
find [13 that the CDW region is the most compressible 
followed by the SF phase, in sharp contrast to a uniform 
CDW which has a gap to charge excitations set by the 
near-neighbor repulsion Vj. Experiments can measure 
5(fc) and therefore r2(fc) ^1] which would serve to verify 
the presence of solitons or other kinds of order. Similarly, 
as commented earlier, the sound velocity is given by the 
linear slope at small fc. 
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